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ABSTRACT 

Elliptical galaxies are modelled with a a 4-component model: Sersic stars, ACDM dark 
matter (DM), a /3-model for the hot gas and a central black hole, with the aim of establishing 
how accurately can one measure the total mass within their virial radii. 

DM is negligible in the inner regions, which are dominated by stars and the central black 
hole. This prevents any kinematical estimate (using a Jeans analysis) of the inner slope of the 
DM density profile. The gas fraction rises, but the baryon fraction decreases with radius, at 
least out to 10 effective radii (i? e ). Even with line-of-sight velocity dispersion (VD) measure- 
ments at 4 or 5 R e with 20 km s _1 accuracy and perfectly known velocity anisotropy, the total 
mass within the virial radius (r v = r2oo) is uncertain by a factor over 3. The DM distributions 
found in ACDM simulations appear inconsistent with the low VDs measured by Romanowsky 
et al. (2003) of planetary nebulae between 2 and 5 R e , which typically (but not always) imply 
very low M / Ls, and a baryon fraction within r v that is greater than the universal value. Re- 
placing the NFW DM model by the new model of Navarro et al. (2004) decreases slightly the 
VD at a given radius. So, given the observed VD measured at 5 R e , the inferred M/L within 
r v is 40% larger than predicted with the NFW model. Folding in the slight (strong) radial 
anisotropy found in ACDM (merger) simulations, which is well modelled (much better than 
with the Osipkov-Merritt formula) with (3(r) = 1/2 r j (r + a), the inferred M/L within r v 
is 1.6 (2.4) times higher than for the isotropic NFW model. Thus, the DM model and radial 
anisotropy can partly explain the low PN VDs, but not in full. The logarithmic slope of the 
VD at radii of 1 to 5 R e , which is insensitive to radius, is another measure of the DM mass 
within the virial radius, but it is similarly affected by the a priori unknown DM mass profile 
and stellar velocity anisotropy. 

Some of the orbital solutions produced by Romanowsky et al. (2003) indicate that 
NGC 3379 has a dark matter content at least as large as cosmologically predicted, and the 
lower Mj Ls of most of their solutions lead to a baryonic fraction within r v that is larger than 
the universal value. In an appendix, single integral expressions are derived for the VDs in 
terms of general radial profiles for the tracer density and total mass, for various anisotropic 
models (general constant anisotropy, radial, Osipkov-Merritt, and the model above). 

Key words: galaxies: elliptical, lenticular, cD — galaxies: haloes — galaxies: structure — 
galaxies: kinematics and dynamics — methods: analytical 



1 INTRODUCTION 

Whereas much work has been devoted to constraining the distri- 
bution of dark matter in spiral galaxies from a multi-component 
(disk, bulge, halo, and gas) modelling of their rotation curves (e.g. 
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IPersic et al]|l99rjt ISalucci & Burkertll200fl) . there have not been 
analogous analyses for elliptical galaxies, where the components 
are stars, dark matter, hot gas and the central black hole. Such an 
analysis is much more difficult for elliptical galaxies, which con- 
trary to spirals, have little rotation llllinsworth 1977), so that one 
cannot directly infer the mass profiles from circular velocities (as- 
suming nearly spherical mass distributions). 

Instead, one has to analyse the velocity dispersions as a func- 
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tion of position, and this analysis involves solving the Jeans equa- 
tion, which in spherical symmetry, assuming no streaming motions 
(including rotation), is 



dr r r 



(1) 



where £ is the luminosity density of the galaxy, a r its radial velocity 
dispersion, and where the anisotropy parameter is 



(2) 



with at — <Je = &<j> the ID tangential velocity dispersion, so that 
/3 = corresponds to isotropy, /3 = 1 is fully radial anisotropy, 
and /3 — * — oo is fully tangential anisotropy. 

The luminosity density is easily obtained by deproject- 
ing the surface brightness profile, and the situation is sim- 
plified by the recent consens us on the applicabilit y to virtu- 
ally all elliptical galaxies ICaon, Capaccioli. & D'Onofrioin"993t 
iBertin. Ciotti. & Del Principe! 1200 2 11 of the gene ralisation (here- 
aft er, Sersic law ) of the R 1 ^ law Ide V aucouleur sll948h proposed 
bv lSersid fl968). which can be written as 



I(R) = 7o exp 



(R_\ 1/r 
\as ) 



(3) 



where / is the surface brightness, as the Sersic scale parameter, 
and m the Sersic shape parameter, with m = 4 recovering the 7? 1 / 4 
law. Moreover, strong correlations have been reported between the 
shape parameter m an d either luminosity or effective (ha lf pro- 
jected light) radius R e JCaon et alll993HPrugniel & SimierJl997l 
and references therein : lGraha^n^^Guzm^rl2O03T) !'^ 

A serious difficulty in modelling the internal kinematics of el- 
liptical galaxies and clusters of galaxies, considered to be spherical 
non-rotating systems, is that the Jeans equation Q has two un- 
knowns: the radial profiles of total mass distribution and velocity 
anisotropy and this mass / anisotropy degeneracy requires some as- 
sumption on anisotropy to recover the total mass distribution. One 
solution is to analyse the veloci ty profiles or at least th e 4th or- 
der velocity moments fMerrittll 19871: iRix & WMtdll992t [Gerhard 
arl2 



IE,. 

ll993tlLokas & Mamorl2003tlKatgert. Biviano. & Mazurel2004 

One can go further in modelling the internal kinematics of 
elliptical galaxies or clusters of galaxies with the Schwarzschild 
1 1979) orbit modelling method, which supposes a form of the 
potential and minimises the differences between the observations 
(maps of surface brightness, mean velocity, and velocity dispersion, 
or even of velocity profiles) and their predictions obtained by suit- 
able projections of linear combinations of a set of orbits that form 
a basis in energy - angular momentum phase space. One can then 
iterate over the form of the potential to see if one finds significantly 
better fits to the observational data. Similarlv. lGerhard et alJ ll998f) 
suggested working with a set of distribution functions, which one 
hopes forms an adequate basis, and one then minimises the linear 
combination of the distribution functions that match best the ob- 
servations. However, such recent analyses rarely place useful con- 
straints on the gravitational potential (with no priors), with the ex- 
ception of Kronawitter et al. (2000), who find that some elliptical 
galaxies are consistent with constant mass-to-light ratios while oth- 
ers show mass rising faster than light. 

The key issue is to obtain more di stant tracers of the g ravi- 
tational potential. Two su ch studies, by Mendez et alJ 1200 it) and 
iRomanowskv et all l2003h . of the distribution of the line-of-sight 
velocities of planetary nebulae (hereafter, PNe) on the outskirts of 
a total of four moderately luminous (L*) nearby elliptical galax- 



ies indicate fairly rapidly decreasing PN velocity dispersion pro- 
files. IRomanowskv et al] carefully analyse one of their galaxies 
(NGC 3379) with the Schwarzschild method, and their favourite 
conclusion is that the dark matter content of these galaxies is low at 
5 effective radii, and extrapolates to a very low mass-to-blue-light 
ratio of T_g = 33 at 120 kpc, which they consider as the virial 
radius. This low mass-to-light ratio at the virial radius (hereafter, 
virial mass-to-light ratio) appears to be at strong variance with the 
cosmological predictions: Two recent works suggest that Tb has a 
non-monotonic variation with mass (or luminosity), with a min- 
imum value arou nd 80-100 for luminosities Lb ~ 10 9-11 Lq 
jMarinoni & Hudson! 120021 : lYang et alJ 120031) . both studies pre- 



x 10 10 L 



dietin g Tb ~ 100 /i7 o for Lb = Lb,* = 
jMamon & Lokasl2005l) . 

There are various alternatives to estimating the mass profiles 
of elliptical galaxies through the Jeans equations, in particular mod- 
elling the X-ray emission arising from hot gas assumed in hy- 
drostatic equilibrium in the gravitational potential, through the ef- 
fects of gravitational lensing and through the kinematic analysis of 
galaxy satellites. 

The analyses of diffuse X-ray emission in elliptical galax- 
ies have the advantage that the equation of hydrostatic equilib- 
rium, which is the gas equivalent of the Jeans equation Q, has 
no anisotropy term within it, so in spherical symmetry, one can 
easily derive the total mass distribution. However, the derivation 
of the total mass profile requires measuring the temperature profile 
and its gradient, and unfortunately, even with the two new gen- 
eration X-ray telescopes XMM-Newton and Chandra, it is diffi- 
cult to achieve such measurements beyond half the virial radius for 
galaxy clusters IPointecouteau et alll2005l : IVikhlinin et alj|2005l) . 
and even much less for elliptical galaxies. Moreover, the X-ray 
emission from elliptical galaxies is the combination of two com- 
ponents: diffuse hot gas swimming in the gravitational potential 
as well as direct emission from individual stars, and it is highly 
difficult to disentangle the two (see lBrown & Bregmanll20(3lT) . As- 
suming tha t all the X-ray em ission is due to the diffuse hot gas, 
lLoewenstein & White! 1 19991) provide interesting constraints on the 
dark matter in luminous (L > L*) ellipticals, with a dark matter 
mass fraction of > 20% (39%) at R = R e (5 Re). 

Weak gravitational lensing is yet another avenue to analyse 
the gravitational potential of elliptical galaxies. Since the signal 
is much too weak for individual galaxie s, one has to stack thou- 
sands of galaxies together. In this manner, IWilson et alJfcOOll) find 
that the gravitational lensing shear falls off with angular distance 
as expected for a structure whose circular velocity is roughly in- 
dependent of radius, which suggests appreciable amounts of dark 
ma tter at large radii . More over, from their analysis of SDSS galax- 
ies. iGuzik & SeliaklfcOOA find that L, galaxies have mass-to-light 
ratio T of order of 100 in the B band, again implying substan- 
tial dark matter, since the stellar contribution to the mass-to-light 
ratio is thought to be only Y * ~ 8 in the blue band (see e.g . 
iGerhard et aUbOOll) . Recently. iHoekstra. Yee. & Gladdersl 120041) 
derived from the gravitational lensing shear profile of RCS galax- 
ies, a lNavarro. Frenk. & White! (hereafter. NFW) mass within the 
virial radius, corresponding to Tb ~ 59 ± 5. 

The constrai nts from th e internal kinematics of galaxy satel- 
lites, pioneered bv lZaritskv & WhiteNl994l). are now showing con- 
sistency with the ACDM models iPrada et all2003l) . However, one 
needs to stack the data from many galaxies, so errors in stacking 
can accumulate, and the method is very sensitive to the correct re- 
moval of interlopers IPrada etHll . 

The constraints on the dark matter are greatly strengthened 
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when combining the internal kinematics with either the X-ray 
or gravitational l ensing approach. Combining internal kinemat- 
ics with X-rays, lLoewenstein & White! 1 19991) are able to con- 
strain the dark matter contents of ellipticals, alth ough they make 
use of the n ot very realistic Osipkov-Merritt I OsipkovJ 1 19791 : 
lMerritJfl985l) anisotropy (see Fig. |2| below). Combining internal 
kinematics with the cons t raints from strong gravitational lensing, 
ITreu & KoopmansI l20oil2004l) are able to constrain the slope of 
the inner density profile of the dark matter. 

On the theoretical side, large-scale dissipationless cosmo- 
logical A-body simulations have recently reached enough mass 
and spatial resolution that there appears to be a convergence 
on the structure and internal kinematics of the bound struc- 
tures, usually referred to as halos, in the simulations. In par- 
ticular, the density profiles ha ve an outer slope of ~ —3 
and a n inner slop e between —1 jNa varro. Frenk, & White 
I l996l) and -3/2 jFukushiee & Makinolll997l: iMoore et al 
Ghigna et al. 2000). In this paper, we consider the general formula 
that ljing & Suto 1 2000) found to provide a good fit to simulated 
halos: 



1995 



1999; 



p(r) oc ( — ) 



1 + 



(-) 

\a.dJ 



(4) 



where the absolute value of the inner slope a = 1 or 3/2, and 'cf 
stands for dark. These profiles fit well the density profiles of cos- 
mological simulations out to the virial radius r v , wherein the mean 
density is A « 200 times the critical density of the Universe, and 
are characterized by their concentration parameter 



0,1 



(5) 



Recently, a number of numerical studi es have prop osed better 
analytic fits to the radial pr ofiles of den sity ( Die mand et all2 004a). 
density logarithmic slope l Navarro et al. 2004) or circular veloc- 
ity |Sto^etail200l|sto^200^ profiles of simulated halos. In 
particular, the formula of Navarro et alJ is attractive because it con- 
verges to a finite central density at very small scales and has an 
increasing outer slope with a finite mass. Moreover, fitting to the 
logarithmic slope of the density profile is to be preferred on the 
grounds that fitting the density or circular velocity involves fitting a 
single or double integral of the density slope, thus possibly missing 
details smoothed over in the integrals. We shall therefore focus on 
the density profile o flNavarro et al] (hereafter. Nav04) in this work. 

Two recent studies have placed interesting constraints 
on the dark matter distribution within elliptical galaxies. 
iBorriello. Salucci. & Danesel 120031) compute the aperture velocity 
dispersions of models with an NFW dark matter component and a 
Sersic luminosity component. They find that the dark matter con- 
tent must be either low or with a concentration parameter roughly 
3 times lower than fit in the dissipationless cosmological TV-body 
simulations, for otherwise the fundamental plane of elliptical galax- 
ies is curved more than allowed by observations. Na politano et alJ 
(2005) estimate the mass-to-light ratio gradient for a large number 
of observed ellipticals. They find that very luminous ellipticals ap- 
pear to have large M / L gradients consistent with important quanti- 
ties of dark matter within the virial radius, while the lower luminos- 
ity L < L* ellipticals appear to have very low gradients, that can 



1 For the standard ACDM parameters, Q m = 0.3, Q\ = 0.7, one finds 
lEke et alJll99dlLokas & Hoffmanll200ll) A = 102, but many cosmolo- 
gists prefer to work with the value of 200, which is close to the value of 178 
originally derived for the Einstein de Sitter Universe (Q m = 1, Q\ = 0). 



be explained either by a strong lack of dark matter, or else by a dark 
matter concentration parameter much lower (c ~ 1 — 6) than ex- 
pected from cosmological TV-body simulations (c ~ 15). It re mains 
to be seen how dependent are the conclusions of iNapolitano et all 
on the kinematical modelling of their ellipticals. 

Our basic goal is to place constraints on the total 
mass d istribution of elliptical galaxies. In a companion paper 
I Mamon & Lokas 2005, hereafter, paper I), we show that the NFW, 
Jine & Suto (hereafter, JS-1.5), and Nav04 density profiles found 
in dissipationless cosmological simulations cannot represent the to- 
tal matter distribution in elliptical galaxies, because: 1) the local 
mass-to-light ratios are then far below the generally admitted val- 
ues for the stellar component; and 2) the aperture and slit-averaged 
velocity dispersions are much lower than observed for their lumi- 
nosities (i.e. through the lFaber & Jacksonll 19761 relation). We also 
show that the highly concentrated NFW models that X-ray ob- 
servers have fit to the total mass density profile of ellipticals are 
an artifact of fitting to the combination of a Sersic stellar compo- 
nent and an NFW dark matter component, and that the fits cannot 
be very good. Finally, we argue that the very low local mass-to-light 
ratios and aperture velocity dispersions found imply that the stellar 
component fully dominates the internal kinematics of the inner few 
effective radii, suggesting that there is little hope in recovering the 
inner slope of the dark matter density profile. All these conclusions 
suppose that the total mass density profile inside an effective radius 
does not sharply steepen to a slope of order —2 or steeper. 

In the present paper, making no direct assumptions on the total 
mass density profile, we focus on the velocity dispersions and mass 
profiles in the outer parts of elliptical galaxies, from 1 to 5 effective 
radii. For this, we build a detailed 4-component model of elliptical 
galaxies, with stars, dark matter, hot diffuse gas and a central black 
hole (as we don't know, a priori, whether these last two compo- 
nents affect the kinematical analysis of dark matter in ellipticals). 
It is very possible that the dark matter reacts to the presence of 
the baryon component, which can be treated by the approximation 
of adiabatic contraction, and recent simulations with gas suggest 
that the resulting dark matt er density p rofile scales as 1/r 2 over a 
wide range of radii l Gnedin et al. 2004 using cosmological simu- 
lations, and lDekel et alj|2005l using galaxy merging simulations). 
Whereas such an "isothermal" distribution is a real possibility for 
the dark matter, it remains u nclear how fee dback processes from 
th e baryons (as modelled by ICox et alj|2005l for the results given 
bv lDekel et alJ) may affect this result. We therefore assume here, 
as in paper I, that the dark matter distribution follows the more ro- 
bust predictions from dissipationless ACDM cosmological A-body 
simulations. 

The plan of our paper in as follows: In Sec.|2| we describe 
our 4-component model, in Sec. 13. II we consi der the mass-t o-light 
ratios consistent with the orbit modelling of Roma nowskv et all 
(2003) on the stellar and PN kinematics in NGC 3379, and in 
Sec. 13.21 we estimate the velocity anisotropy of the particles in 
dissipationless cosmological A-body simulations. We study the ef- 
fects of velocity anisotropy in Sec. 14.11 compare the importance 
of the different components in Sec. 14.21 and ask whether we can 
weigh the dark matter component in Sec. 14.31 We discuss our re- 
sults in Sect. [5] and r eflect on the ve ry low mass-to-light ratio re- 
ported by Romanows kv et alJ 12003) for intermediate luminosity 
elliptical galaxies. In Appendix lAl we derive single quadrature ex- 
pressions for the radial profiles of the line of sight velocity disper- 
sion in terms of general tracer density and total mass profiles for 
four simple anisotropy models. 
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2 BASIC EQUATIONS 

We now highlight the 4-component model that we adopt for el- 
liptical galaxies: stars, diffuse dark matter, hot gas, and a central 
black hole. We neglect any contribution from dark baryons (e.g. 
MACHOs). 



2.1 Distribution of optical light 

We begin with the distribution of optical light, referring the reader 
to paper I for details. The Sersic (eq. (3), see detailed properties 
in iGraham & Dri veil 12005) optical surface brightness profile that 
represents the projected stellar distribution, can be deprojected ac- 
cording to the approximation first proposed by Prugniel & Simien 
( 1997) 



£(r) 
t{x) 

li 



lil(r/as) 
axp ( 
T(2m) 



— v ( l/m\ 

x exp ( — x 1 ) 



(6) 
(7) 

(8) 

(9) 



T[(3-p)m]j 2 a s ' 
p ~ 1.0 - 0.6097/m + 0.05463/m 2 , 

where the last equation is from L ima Neto et alj fl999). 

The integrated luminosity corresponding to equations J6j, 
and J8} is then iLima Neto et alJ) 



7 [(3-p)m,x 1/m ] 



L 3 (r) = LL 3 (r/a s ) 

r[(3-p)m] 
where the total luminosity of the galaxy is 

L = 2-Km r(2m) I a% = iirm Y [(3 - p) m] ix a% 



(10) 
(11) 

(12) 

as obtained by Youn g & Curriel < 19941) from the Sersic surface 
brightness profile of equation {3J, and which matches exactly the 
total luminosity obtained by integration of Lima Neto et al.'s ap- 
proximate deprojected profile. 

It is useful to express radii in terms of the effective radius, R e , 
which is the radius of half projected light, where 

Re = b m a s (13) 
b ~ 2m- 1/3 + 0.009876/m , (14) 

where the latter relation is from lPrugniel & Simienl 1 19971) . In pa- 
per I, we showed that R e and m are fairly well correlated with total 
luminosity: 



log h 7Q R, 



(fit) 



0.34 + 0.54 log Lio + 0.25 (log Lio) , (15) 



logm (flt) = 0.43 + 0.26 logL 10 - 0.044 (logL 10 ) 2 , (16) 

where Lio = h% Lb/ (l0 10 L©), Re is measured in kpc, and 
with Ho = 100 h = 70 hro kms' 1 Mpc" 1 . 2 Then, equations O 
and (14} lead to 



as 



R, 



(fit) 



(17) 



2.2 The central black hole 



High spatial resolution spectroscopic studies of ellipticals have 
shown that they almost always harbour a supermassive black hole 

2 Unless noted otherwise, we adopt h,7o = 1. 



of ma ss 0.2% jFaber et all 19971) . 0.3% jKormendv et alll997l) . or 
0.6% lMagomariTHdT998l) that of the stellar component. We thus 
define the fraction of black hole mass to stellar mass 



9BH 



T*L 3 (rg 



(18) 



We could have defined gBH in terms of the total luminosity in- 
stead of that at the virial radius, but the two differ by typi- 
cally less than 0.1%, which is much less than the uncertainty 
on 3bh, and the latter s caling generat es simpler equations be- 
low. A recent analysis I Haring & Rix| 120041) favours <?bh = 
0.15% (Mb u ig C /6 x 10 11 Mq)°' 12 , and we will adopt below 
g Bn — 0.15%. The precise value of <?bh has a negligible effect 
on the constraints on dark matter in elliptical galaxies. 



2.3 Scalings of global properties 

We adopt a fid ucial luminosit y of L = L». In paper I, we de- 
rived, from iLiske et alj J2003I) . a blue-band luminosity L,.b = 



1.88 10 1 



'■70 



Lq, i.e. an abso lute magnitude of Mb = —20.24, 



using (e.g.. lColina etalll996l) 



Mf = 5.45 . 



(19) 



Our choice of Lb translates to m = 3.15 (eq. 1161 ). R e = 
3.2 kpc (eq. 1151 ) and as = 11.6 h^pc (eq. 1171 ). 

In paper I, we argued that the blue stellar mass-to-light ratios 
of elliptical galaxies lay in the rough range from 5 to 8, where the 
uncertainty is caused by the uncertain initial mass function and its 
lower and upper mass cutoffs, as well as metallicity and stellar evo- 
lution code. In the present paper, we adopt the mean: T»,s = 6.5, 
unless specified otherwise. 

We also showed in paper I that the mass-to-light ratio of the 
Universe is 



Tb — f^m T closure, B — -z — . ° — 390 /170 

8nGjB 



(20) 



given the luminosity density found by Liske et al. I2003J) and 
Blanton et al. 1 2003) in the 2dFGRS and SDSS surveys, respec- 
tively. 

We define the mass-to-light ratio bias 



h T 
T 



(21) 



where T = Mt /L(r v ), with Mt the total mass within the virial ra- 
dius. If the Universe is unbiased, i.e. the mass-to-light ratio within 
the virial radius will be Tb = Tb = 390 hro - The internal kine- 
matics of galaxy clusters are consistent with the universal mass-to- 
light ratio (e.g. lLokas & MamonlEooH deriv e T B = 351 for the 
Coma clust er). Followi n g the predictions of M arinoni & H udson 
(2002) and lYang et al] J2003I) . we will adopt a standard value 
T B = T s i d = 100/170. 

We also define the baryon fraction bias 



b b =t = 



fb 



f (Qb/Qm) 



(22) 



where fb is the baryon fraction within the virial radius, while 
f b = Qb/^m — 0.14 is the mean baryon fraction of the Uni- 
verse. We used here the bi g-bang nu cleosynthesis measurement 
Q b = 0.02 h~ 2 = 0.041 ]o'MearaetaiTl200ll) whi ch is con- 
sistent with the value obtained by ISpereeletalJ T2003) from the 
WMAP CMB experiment, and fi m = 0.3. 



© 2005 RAS, MNRAS 000, 000-000 



Weighing dark matter in ellipticals 5 



Given the fractions of mass in stars, the central black hole, hot 
gas and dark matter, all at the virial radius: 



/* = 

/bh = 

fd = 

h = 



9BH /* , 

i-f b = i-b b 7„, 

!-(/. + /bh + /a) 



(23) 

(24) 
(25) 
(26) 



(27) 



where we assumed that the central black hole originates from 
baryons (this assumption has a negligible effect on what follows), 
it is easy to show that these two biases are related through 

l +g _+ gB H (T.\1+9(T. 

where g is the gas to star ratio within the virial radius, and T„ is the 
mass-to-light ratio of the stellar population, assumed independent 
of radius (in confor mance with the v ery weak colour gradients in 
ellipticals, e.g. lGoudfrooii et alll994h . 

If mass is biased relative to luminosity in elliptical galaxies 
(br 7^ 1), we can assume that the baryon fraction is unaffected 
(bt = 1, and g must vary according to eq. 1271 ). Indeed, it is dif- 
ficult to conceive of a mechanism that will segregate baryons from 
dark matter within a radius as large as the virial radius. For the gas 
to star ratio to be positive, one then requires 

(¥) • 



T s > =48 



h 

In other words, considering more mass in stars requires a larger 
minimum total mass at the virial radius for the baryon fraction to 
retain its universal value. Turning the argument around, a low total 
mass-to-light ratio at the virial radius would imply that the baryonic 
fraction is greater than the universal value, i.e. bb > 1. Indeed, the 
general equation 1271 implies 

where the equality is for a negligible gas-to-star ratio g. 
2.4 Distribution of dark mass 

We consider here three dark matter models: the NFW model 
with inner s l ope — 1, the generalized NFW model introduced by 
Ijing & Sutol l2000|) with inner slope — 3/2 (JS-1.5), and the con- 
vergent model of iNavarro et alJ 120041 Nav04). The dark matter 
density profile can generally be written (see paper L especially for 
the Nav04 model): 



Pd(r) = -^f-r (y^s) Pd(r/a d ) , 
g d {c) \4:TTri/ 



9d(c) 

y~ a (i + y) c 



(30) 



a = 1 (NFW), a = 3/2 (JS-1.5) . 



Pd{y) = 



exp (-2/iy 1/fl ) (Nav04). 
(]n(y + l)-y/(y + l) 



9d(y) = <2 



sinh 1 y/y - y/y/(y + l) 



(NFW) . 



(JS-1.5) . 



(31) 



(32) 



k | (2/i) 1 - 3 ^ 7 [3 M ,2M2/ 1/A1 ] (Nav04), 

where a is the absolute value of the logarithmic slope of the inner 
density profile (for the NFW and JS-1.5 models), c is the concen- 
tration parameter (eq. Q), ad is the radius where the logarithmic 



slope is equal to -2 (NFW, Nav04) or -9/4 (JS-1.5, for which 
ad/2 is the radius where the slope is — 2), ji ~ 6 (paper I), and 
sinh - x = ln(a; + s/x" 2 + 1) for x > 0. In equation 1301 . Md is 
the dark mass within the virial radius, defined such that the mean 
total density within it is A = 200 times the critical density of 
the Universe, p CT j t = 3 Hq/(8wG), yielding a virial radius (see 
lNavarroetalll997l) 



2 G My 



1/3 



163 h kpc 



= 206 h 7 g kpc 



hM t 
10 12 M Q 

h 70 Mt 



1/3 



1/3 



(33) 



10 12 Mq 

for A = 200. |jing & Sutol (;2000) measured the concentration pa- 
rameter c (eq. Q) from their ACDM (with cosmological density 
parameter Q m = 0.3 and dimensionless cosmological constant 
Qa = 0.7) simulations, which can be fit by the relations 



10.2 M 



4.9 Mi 



12 
0.13 



(NFW) . 
(JS-1.5) 



(34) 



where M 12 = hM t /10 12 M Q . In paper I, we derived the concen- 
tration parameter of the Nav04 model: 

C = 8.1Afr 2 ail -°- 0151ogMl2 (35) 

which is similar to the concentration parameters in equations 1341 . 
For Lb = L, t B, equations I15i and <33> lead to 



Re 



79 (1*) 

Vioo; 



1/3 



(36) 



The cumulative mass of the dark models used here can all be 
written 



M d (r) = 
M d (y) = 



M d M d (r/a d ) 
9d{y) 



(37) 
(38) 



9d{c) ' 

where gd(y) is given in equation <32t . 3 

A complication arises because the luminosity at the virial ra- 
dius is very close to but not exactly equal to the total luminosity. 
We write out the total mass at the virial radius, both as a mean den- 
sity threshold, and as the luminosity at the virial radius times the 
mass-to-light ratio: 



M t = — 4 A p C ri t 



2G 



T L L 3 (r v /a a ) , 



(39) 



(where we used eq. 12 II ) and solve the last equality of equation J39t 
for r v (the correction on r v turns out to be negligible as the lumi- 
nosity has nearly fully converged at the virial radius). 



2.5 Distribution of gas mass 

Elliptical galaxies also shine X-rays. Both discrete sources and a 
hot diffuse interstellar medium contribute to their X-ray emission, 
and the latter component should contribute non-negligibly to the 



3 Note that the definition of gd is the inverse of the definition of g given by 
lLokas & Mamonl 1 200 1 ) for the NFW model. 
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mass budget of ellipticals. Indeed, if elliptical galaxies are unbi- 
ased matter tracers (6t = 1), then according to equation I23t . 
stars contribute at the virial radius to a fraction of the total mass of 
/, = T»/T = 0.017 (T*/6.5). This is much less than the baryon 
fraction at the virial radius assuming no baryon bias (6b = 1), i.e. 



fb 



ft = 0-14 ft™ 



(Sect. |23j- So, if ellipticals have the same 



baryon fraction within their virial radius as the full Universe, then 
the hot gas in ellipticals accounts for 0.14 — 0.025 — 0.11 
(for hro = 1) of the total mass within the virial radius (assum- 
ing no dark baryons), which is 0.11/0.025 ~ 4 times more than 
stars. If elliptical galaxies have virial mass-to-light ratios as low 
as Tb = 100, then the gas-to-star ratio is lower (for fixed bary- 
onic fraction). According to equation J27t . for no baryon bias 
(&6 = 1), the gas dominates the stars at the virial radius (g > 1) for 
T B > 2T«. B /7, = 93 ( T«. b / 8). 

brown & Bregma n (2001) have fit the X-ray surface bright- 
ness profiles of luminous ellipticals (jointly with a component for 
discrete sources following the i? 1//4 law), and find that the usual so 
called /3-model 



Pair) = pop g (r/r c ) (40) 
p„(z) = (l + , 2 )~ 3 ^ /2 , (41) 

provides a good representation of the distribution of the hot gas, 4 
with (P g ) = 0.49 and 

Re 

r c ~ — , 

q 



(42) 



where q ~ 10. Usi ng a larger sample of elliptical galaxies, 
lO' Sullivan et alj l2003h find a mean /3 g of 0.55, but do not esti- 
mate r c /R e . These parameters should be considered tentative, as 
the subtraction of the hard stellar component is quite uncertain. 

As for the other non-stellar components, the gas enters our 
analysis only through its cumulative mass (see the Jeans equa- 
tion Q). Adopting 

Pa = I/ 2 . (43) 
the cumulative mass distribution arising from equation I40i is, 



M g {r) = 4tt p rlM g {r/r c ) 



M B {z) = -z 



(1 + 2 



2x1/4 



(1 3 3 
V2' 4' 2' 



(44) 
(45) 



where 2 F 1 (a, b; c; as) = r(c)/[T(6)r(l-A)] £ t^jl-t)^ 1 (1- 

tx) ~ a is a the Hypergeometric function. Note that M g (l) ~ 0.239. 
To an accuracy of 2.7% for all z, one has 

-l/ 7 



M g (z) 



(46) 



with 7 = 2 1/8 ~ 1.0905. 

We write the gas mass profile as 



M g {r) = M g {r v ) 



M g (r/r c ) 
M g (r v /r c ) 



(47) 



and normalise the gas component with the baryon fraction within 
the virial radius: 

M (r 1 _ T 

-^f 1 = fb ~ (1 + ,9bh) U = h f b - (1 + 3 bh) ^ . (48) 



4 Our parameters z and f3 g have nothing to do with a redshift or velocity 
anisotropy, respectively! 



The divergence of the gas mass profile at large radii is very 
severe (M g oc r 3 ^ 2 ), leading to a divergent gravitational potential, 
and to a local gas fraction much too large. We therefore assume 
that, beyond the virial radius, the ratio of local baryon (gas + stellar) 
mass density to local total matter density is equal to the universal 
baryon ratio f b — f2(,/O m , yielding 



Pai r ) 



h 



fx 



■p d {r) -TUir) 



for r > r v , 



(49) 



where pd is the dark matter density. 

The local gas density profile will therefore be discontinuous 
at the virial radius. Physically, one expects a shock to occur at the 
interface between the gas infalling into the galaxy and the gas in 
equilibrium within the galaxy, and this shock should occur very 
close to the virial radius. With equations j40}, E), {42), j^), <44l . 
<46t . and 1481 . the local gas density profile just within the virial 
radius is (dropping the negligible 1 + qbh term in eq. 1481 ) 



Pair-o ) 



(hf. 



- T» 



[l + rl/rl) 



-3/4 



( 



\ 47rrj 



(50) 



M g (r v /r c ) 

while with equations 1301 . 13 li . 1321 . and I49i . the local gas density 
just outside the virial radius is 



Pair v ) = 



[c/(c+i)r 



' ft, 



1 - h f b 



1-/6 



(r„/as) 3 tirv/as) 



T* 



( 



M, 



(51) 



mr[(3- P ML 3 (r„/a s ) TJ \4irri 

According to equations I50I and <51> . the discontinuity of the 
gas density at the virial radius can be written 



Pgjrj) 
Pai r v) 



[1 + zl) 3/4 M g {z v ) 



b b Tf b -T, 



T/ 6 



l-hf h \ [c/(c+l)] J 
l-7i 



9d(c) 



mr[(3-p)m] L 3 {x v ) 



(52) 



wherea;^ = r v /as andz v = r v /r c . For a wide variety of plausible 
parameters T„, &(,, T, a, c, x v and z v , equation <52t yields density 
ratios p 9 (f^")/p s (r^") between 2 and 3.5, in accordance with the 
standard Rankine-Hugoniot conditions (density ratio smaller than 
4). 

The cumulative gas mass beyond the virial radius is then 



M g (r) = M g (r v ) + 



fb 



-T. [L(r 
yielding a normalized gas mass satisfying 



>— [M d (r) - M d {r v )\ 
Hr v )] , 



(53) 



M g jz) 
M g (z v ) 



M g jr) 
M g ir v ) 

7 b 



k 

fa 



(1 

f* 
fa 



-ft) 



M d ( JL) 
\a d J 



- 1 



Ls jr/as) 
La, ir v /a s ) 



(54) 



where the second term is negligible (as the integrated luminosity 
has almost fully converged at the virial radius). 
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3 OTHER ASTRONOMICAL INPUTS 



3.1 Lower bound on the virial mass-to-light ratio 

Our lower bound on the mass-to-lig ht ratio within the virial ra- 
dius is taken from the work of Romano wskv et all 1200 3l) . who per- 
formed, for the nearby giant elliptical, NGC 3379, an orbit mod- 
elling of the velocity field of the planetary nebulae combined with 
the radial profiles of the stellar surface brightness and line-of-sight 
velocity dispersion. Although they quote a mass-to-light ratio at 
1 20 kpc, which they estimate as the virial radius, of Tb = 33 ± 14, 
we find 



A(r) = 2l|^=507 : 



(55) 



for r = 120 kpc, T B = 33, and L B = 1.51 x 10 10 L© (from 
iRomanowskv et alT s M B = -20.0 and eq. (T5)). So clearly the 
radius of overdensity 200 is much larger than 120 kpc, and since 
dark matter is believed to be more extended than luminous matter 
(e.g., Fig.|4j, the mass-to-light ratio should be larger than 33. 

We extrapolate lRomanowskv et alt s mass-to-light ratio out to 
the virial radius, as follows. In their orbit modelling of NGC 3379, 
IRomanowskv et all (supporting online material) w rite the density 
profile as the sum of a luminous iHernauisl ll990l) component and 
a dark NFW component: 



p(r 



+ 



2nGr(r + a) 3 4TrGr(r + r 



(56) 



Writing r = ax and r s = ax„, it is easy to show that the mass 
enclosed within radius r satisfies 



GM (r) 



■(57) 



Expressing the mean density at the virial radius as A times 
the critical density of the Universe, one needs to solve for x v — 
r v /a = cx s the equation 



XX-u + U \V*J 



In 1 + 



Xd ~\- X s 



(58) 



where r v i r is the virial radius for the total matter. 

Figure Q shows the resulting values of the concentration pa- 
rameter of the NFW dark matter model versus the mass-to-light ra- 
tio evaluate d at the virial radii, r 2oo (open circles) and rio2 (filled 
circles), for IRomanowskv et alt s 15 orbit solutions (with the val- 
ues of r s , v», and vl/vl given in their supporting online mate- 
rial). Here, we first solved equation 1581 . and substituted the di- 
mensionless virial radius, x v for x in equation <57> , We also plot 
(as diamonds) in Fig.Qanother acceptable solution, which A. Ro- 
manowsky kindly communicated to us, leading to Ts as high as 
164 at r2oo a nd 196 at rica. 

Whereas Roman owskv et all quote Ts = 33 ± 14 at 120 kpc, 
we find, using the same input values of r s , v t , uf/v,, assuming 
R e = 36", 5 (from which we infer a Hernquist scale radius of a = 



5 The effective radius adopted by Romanowskv et al. 1 2003 ), origi- 
nates from aperture photometry i de Vaucouleurs et alj Il99ll) . Determi- 
nations from fits of the differential surface bright ness profile over a 
very l arge radial range yield R e = 53" (dg Vaucouleurs & Capacciolj 
M^9l_ Capaccioli et al. 1990), so the effective radius adopted by 
iR^mmowskyTt^lT mav be underestimated by a factor 1.3. 
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Figure 1. Dark matter (NFW) concentration versus total mass-to-light ratio 
at the virial radius, T200 (open circles) and rirj2 (filled circles), as inferred 
from the orbit solutions of Romanowsky et al. (2003, supporting online ma- 
terial, circles) and Romanowsky (2005, private communication, diamonds), 
using equation 1581 . The oblique horizontal line indicates the concentration- 
mass that we infer in paper I from the dissipationless ACDM simulations of 
Navarro et al. (2004), while the vertical line gives the mass-to-light ratio at 
the virial radius from the cosmoslogical prediction of Marinoni & Hudson 
(2002). The hashed regions give an indication of the uncertainties on the 
concentration-mass relation (factor 1.4) and the mass-to-light ratio at the 
virial radius derived by Marinoni & Hudson (factor 1.6). The arrow indi- 
cates where the solutions lead to a baryonic fraction that is higher than the 
universal value. The filled symbols ought to lie within the hatched regions 
and to the right of the arrow. 



1.0 kpc), Ts(120 kpc) to be in the range 19 — 48, with a mean of 
27 and a median of 22. Extrapolating to the virial radius, defined in 
this paper as r2oo, we find T b in the range 20 to 70, with a mean 
of 32 and a median of 25, while for the ACDM virial radius, rio2, 
we get T b between 22 and 82, with a mean of 37 and a median of 
28. 

Figure Q shows that one of the 15 orbit solutions of 
\Romanowskv et al\ has a combination of virial mass-to-light ra- 
tio and concentration parameter in line with the predictions from 
dissipationless cosmological simulations (T_g,io2 = 70 and c 
10). where the c o ncentr ation parameter predictions are those that 
iNapolitano et alJ 120051) rederived from the dis s ipation less cos- 
mological simulations analyzed by Bullock et al. 1 2001 ). None of 
the orbit solutio ns lead to t he ve ry lo w concentration p arameters 
favoured by B ordello et alj J2OO3I) and lNapolitano et all 

We also note that the Hernquist model contribution to the 
mass-to-light ratio at the virial radius is always near T„,s = 5.6. 
According to equation 1281 . this amounts to a minimum mass-to- 
light ratio of Tb = 41 for the baryonic fraction at the virial radius 
to be equal to the univ ersal va lue. Only three of the fifteen orbit so- 
lutions of Roma nowskv et alJ can be modelled with no baryon bias, 
but for their mean value of Tb — 33, which we adopt as our lower 
limit, the baryon bias is (eq. 1291 ) bt = 1.25 with T„,b = 5.6 or 
bb — 1.45 with Y^b = 6.5. Therefore, in syst ems with total mass- 
to-light ratios as low as the mean derived by Roman owskv et a~I\ 
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the baryon fraction within the virial radius is greater than the uni- 
versal v alue, unless the stellar m ass-to-light ratio is overestimated 
(see also N apolitano et all2005l) . One is then left with the difficulty 
of explaining how baryons can segregate from dark matter as far out 
as the virial radius. 



3.2 Velocity anisotropy 

Given the mass / anisotropy degeneracy mentioned in Sect. Q 
it is difficult to estimate the radial variation of the anisotropy 
parameter, (3 (defined in eq. The simplest solution is to 

assume isotropy throughou t the galaxy. In fac t, f or clusters of 
galaxi es, lMemttNl987t) . lLokas & Mamonl fiw%) and lKatgert et alJ 
2004) each derived isotropic orbits throughout, Merritg and 



Katge rtet alJ by examining the global velocity distribution and 
Lokas & Mamon by performing a joint local analysis of the radial 
profiles of the line-of-si ght velocity dispersi on and kurtosis. With 
their Schwarzschild-like Gerhard et al. 1 1998) modelling of ellipti- 
cal galaxies, Kronawitter et al. (2000) and Saglia et al. (2000) find 
(3 < 0.3 within R e /2. 

It is useful to compare these still very uncertain anisotropy 
estimates with those measured in cosmological TV-body simula- 
tions. Figure |2| shows the anisotropy profiles for five recen t analy- 
ses of dissipationl e ss cosmological simulations bvlDiaf erio|{l99'9l). 
Colin et alJj200fj|).|Rasia et alJ<2004l) . lDiemand et alj<2004bl) . and 
Woitak et alJ 120051) . with the latter four converted from r /rioo to 
r/r2oo assuming 7-100/V200 = 1.37, as expected for Nav04 mod- 
els with the masses (at the level of g alaxy clusters) a s simulated. 
Note that the orbital anisotropies of lDiaferid Jl999l) . IColm et alJ 
(2000), and lWoitaketalJ jl005) include the mean radial an d tan- 
gential streaming motions, while that of iRasia et alJ 120041) does 
not. 

The Figure also displays analytical repre sentati ons of the data: 
the model (dotted curve ) fhat lCarlberg etalJ<1997l) fit to the kine- 
matics of CNOC clusters: 



P(r) = 2f3„ 



(59) 



with p m = 0.65 and r a = 2 r v , which lColfn et alJ 120001) found to 
fit well the anisotropy of the subhalos in their simulation, as well 
as an anisotropy model that also appears to fit well the simulation 
data: 



1 r 
B(r) = - — — 

H( ' 2 r + Ta 



(60) 



for r a = 0.18 (lower solid curve). 

The anisotropy profiles plotted in Figure|2|are for all the parti- 
cles of a halo, instead of for the subhalos of a halo, which are much 
closer to isotropy. Indeed, we are wary of using the subhalos in 
the cosmological simulations, because the number density profile 
of subhalos within halos is suspic ious, as it has a much shallower 
inner slope than the dark matter ( Colin et alJ ll999l) with a nearly 
homogeneous core ( Diemand et al. 2004a), in contrast with the dis- 
tribution of galaxies in clusters ( Carlberg et al. 1997). Interestingly, 
stars show more radial o rbits than dark m atter particles in hydro- 
dynam ical cosmological ISaiz et all2004l) and merger ( Dek el et alJ 
2005) simulations, where the gas is allowed to cool and form stars, 
and, for the latter study, have a feedback effect on the cooling of the 
remaining gas. Dekel et al. find a stellar anisotropy profile that re- 
sembles the model of equation 1601 . with r a — 1.4 R e ~ 0.018 r v 
(using eq. 1361 for the latter approximation). This is shown as the 
higher of the two solid curves in Figure|2| 

Figure [2] indicates that the anisotropy of equation I60i with 



0.6 



0.4 



02. 




0.2 - 



Figure 2. Radial variation of the velocity anisotropy found in cosmologi- 
cal TV-body simulations. The asterisks, triangles, stars, circles, and squares 
correspond to particles within halos of cosmological ACDM TV-body simu- 
lations of Diaferio (1999), Colin et al. (2000), Rasia et al. (2004), Diemand 
et al. (2004b), and Wojtak et al. (2005), respectively. The dotted curve is the 
profile proposed by Carlberg et al. (1997) (eq. 1591 1 and found by Colin et al. 
to fit well the subhalo anisotropies. The solid curves are the analytical fit of 
equation 1601 . with r„ /r v = 0.18 (providing a good fit to the data from the 
dissipationless cosmological TV-body simulations) and 0.018 (matching the 
anisotropy profile of Dekel et al.'s 2005 merger remnants), going upwards 
(red and blue in the electronic version of the journal). The dashed curves 
are the Osipkov-Merritt anisotropy (eq. 1611 ) with r a /r v = 0.25 and 1, 
going downwards. The upper limit is from the analysis of Kronawitter et al. 
(2000) and Saglia et al. (2000), assuming r v /R e ~ 100 as expected for 
our standard parameters (Lg = 10 10 /170LQ and = 100). 



r a = 0.18 r v fit s very well the anisotro py profiles of IColfn et all 
IRasia et all and iDiemand et all <2004bl) from 2% to 10% of the 
virial radius, which happens to be the region where we shall find 
that the anisotropy may play an important role in modelling the 



mass-to-light ratios of elliptical galaxies. Fi* 


;ure|2|also shows the 


commonlv-used OsiDkov-Merritt lOsiDkovl 


1979; 


Merrit!l[L98a) 



anisotropy model, 

r 2 

= ZTT^ > 



(61) 



for r a /r v — 0.25 and 1. Clearly, the Osipkov-Merritt anisotropy 
is a poor fit to the simulations, as it converges to a too high value 
of unity, and worse, decreases to zero too fast at increasingly low 
radii. 

In what follows, we generally assume that the anisotropy of 
the stellar population is equal to the anisotropy of the dark mat- 
ter particle system, but we will also allo w for the much stronger 
radial anisotropy found by iDekel etUll In appendix [A] we de- 
rive the line-of-sight velocity dispersion profile for four different 
anisotropy profiles: constant anisotropy, the limiting purely radial 
case, Osipkov-Merritt, and the anisotropy model of equation 1601 . 
With these single quadratures, one avoids double integration (first 
integrating the Jeans equation to obtain the radial velocity disper- 
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sion, and then integrating along the line of sight to get the line-of- 
sight velocity dispersion). 



4 RESULTS 

If, as we found in paper I, the mass models found in dissipation- 
less cosmological TV-body simulations are not able to reproduce by 
themselves the rather high central velocity dispersions observed in 
elliptical galaxies, this suggests that the central velocity dispersions 
of ellipticals are dominated by the stellar component and by a cen- 
tral supermassive black hole, as we shall see below. We are now left 
with the question of whether the NFW, JS-1.5 and Nav04 models 
are adequate in describing the diffuse dark matter component (ex- 
cluding the central black hole) of elliptical galaxies. 

Dark matter is expected to become significant in the outer re- 
gions of ellipticals, since the black hole affects the central regions, 
and the influence of the stellar component is usually thought to be 
important in the inner regions, at least within R e . We therefore pre- 
dict the line-of-sight velocity dispersion profiles of ellipticals built 
with 4 components: a Sersic stellar component of constant mass- 
to-light ratio T„, a central black hole, a hot gas component, and 
a diffuse dark matter component described by an NFW, JS-1.5, or 
Nav04 model, according to our parameterisations of Sec.|5| We fo- 
cus our analysis on the following set of parameters: a luminosity 
L B = £»,b = 1.88 x 10 10 h~g L Q (implying R e = 3.2 /i" 1 kpc 
from eq. 1151 . m = 3.15 from eq. 1161 . and as = 12 h? pc from 
eq. 1171 ). a stellar mass-to-light ratio T».s = 6.5 and no baryon 
bias (6;, = 1). 



4.1 The effects of velocity anisotropy 

We first check the effects of anisotropy, by computing the line-of- 
sight velocity dispersions of our four-component model, assuming 
different anisotropy profiles. Figure|5|shows the total line-of-sight 
velocity dispersion (i.e., the quadratic sum of the individual veloc- 
ity dispersions of the four components) for four anisotropy models: 
isotropic, radial, and the model of equation I60i for two choices of 
r a . The line-of-sight velocity dispersion profile of the anisotropy 
model of equation 1601 with our standard value of r a = 0.18 r v 
is virtually indistinguishable from that of the isotropic model: 
the anisotropic model producing 2% lower velocity dispersions at 
5 R e . On the other hand, if the velocity ellipsoid is indeed radially 
anisotropic beyond 7^/50, as suggested by the orbital anisotropics 
of the stellar par ticles in the simulations of Sai z et alJ 120041) and 
iDekel et"al] 120051) . then one can obtain a small but non-negligible 
decrement in line-of-sight velocity dispersion: if r a = 0.018 r v , 
one finds decrements of 4% at 2 R e and 9% at 5 Re . Note that with 
Osipkov-Merritt anisotropy, and r a = r v or even r a = 0.25 r v , 
the line-of-sight velocity dispersions (obtained with eq. IA12I ) are 
indistinguishable from the isotropic ones. 

Hence, with an optimistic precision on velocity dispersion of 
5 kms -1 , the slight radial anisotropy found in cosmological N- 
body simulations produces line-of-sight velocity dispersion profiles 
that are virtually indistinguishable within 9 Re from analogous 
profiles assuming velocity isotropy everywhere, while the stronger 
radial anisotropy found for stars in simulated merger remnants dis- 
plays measurable differences beyond 2 R e . 

In comparison, as illustrated in Figure [5] discarding the dark 
matter component has a much stronger effect on the velocity disper- 
sions at a few effective radii than does the strong radial anisotropy 




0.1 1 10 

Figure 3. Line-of-sight stellar velocity dispersion of 4-component ellipti- 
cal galaxy for different anisotropy models (eqs. IA19I and IA16I ), with a 
Nav04 dark matter model, L B = L„ jS = 1.88 X 10 10 L Q , m = 3.15, 
Re = 3.2 h~Q kpc, T», s = 6.5, T B = 100, g BH = 0.0015 and b b = 
1. The short-dashed, solid, and long-dashed curves represent the isotropic 
(see eq. IA1 II ). mildly and sharply anisotropic (eq. 1601 . with r a /r v = 
0.18 and 0.018), using equation IA9I ). respectively, while the dash-dotted 
curve represents a galaxy with no dark matter (T g = fb T^ d = 14) and 
mild anisotropy (eq. 1601 with r a /r v = 0.18) — these four cases are re- 
spectively shown in red, black, blue and green in the electronic version of 
the journal. 



with r a = 0.018 r v , with dispersion 19% and 33% lower than with 
dark matter (at T B = 100). 



4.2 Relative importance of dark matter, stars, gas and the 
central black hole 

Figure [4] shows the contribution of each component to the lo- 
cal mass density, cumulative mass, line-of-sight (with the mild 
anisotropy of eq. 1601 and r a /r v = 0.18) and aperture (isotropic 6 ) 
velocity dispersions. For our standard parameters, Lb = £*,s = 
1.88 x lO lo /if o 2 L , T»,s = 6.5, T B = 100 (the typical value 
inf erred in the cosmo logical analyses of Marinoni & Hudso j2002l 
and lYang et all2003l i.e. b T = 0.26), and g B n = 0.0015, the stel- 
lar component dominates over the dark matter component out to 
R ~ 2 R e for the local density, 5 Re for the enclosed mass, 3 Re 
for the line-of-sight velocity dispersion, while the stars dominate 
the dark matter everywhere for the aperture velocity dispersion. 
In the inner regions, with <?bh = 0.0015, the central black hole 
is the dominant component for R < 0.015 R e for the enclosed 
mass, 0.006 R e for the aperture velocity dispersion, but only out to 
0.003 Re (i.e., ~ 10 pc) for the line-of-sight velocity dispersion. 

6 The aperture velocity dispersion with the anisotropy of eq. 1601 is diffi- 
cult to express in terms of a single quadrature, but the difference with the 
isotropic case should be very small, since the difference between isotropic 
and mildly anisotropic line-of-sight velocity dispersions is negligible, see 

Fig.m 



© 2005 RAS, MNRAS 000, 000-000 



10 Gary A. Mamon and Ewa L. Lokas 



0.01 0.1 1 10 1 

10 




-2 
10 4 

~p 103 L 
J 10 s r 



^ 10 ZL 



5 i 

S3 

o 0.1 



1 1 iiiiii| — 1 1 m i n i — 1 1 iiiiii| — i-yy 



0.01 
200 




1 1 MifHf 1 1 m i n i — 1 1 iiiiii| — n-Hf 




-- s io- 8 r^* 



H IO" 3 h 



f i i iiiiii| — i i iiii n | — i i iiiiii| — i-H-m- 



© 



100 r 



10 

5! 



0) 

> 

3 



1 - 



200 :- 



' w 60 
E 40 

20 



'8 = 

6 - 
4 



10 100 

rrmj 1 — m 





i o.i e^X^ r& 



total 




* . . BH 



gas 



uL 



IXLUJ I l_U 



0.01 0.1 1 10 

R / R e 



100 



Figure 4. Contribution of stars (short dashes), dark matter (light curves; dashed for NFW, dotted for JS-1.5, and solid for Nav04), hot gas (long dashes), 
and the central black hole (big dots) to the local mass density (upper left), cumulative mass (upper right), local (middle left) and cumulative (middle right) 
mass-to-light ratios, mildly anisotropic (eq. 1601 with r a = 0.18 r v ) line-of-sight velocity dispersion (lower left) and isotropic aperture velocity dispersion 
(lower right); for a giant elliptical with Lg = L t g = 1.88 X 1O 1o /i^"q 2 L0 (corresponding torn = 3.12, R e = 3.2 h^n kpc and ag = 11.6 h™ 1 pc), 
Y* g = 6.5, Tg = 100 (i.e., by g = 0.26), yielding r v /R e = 79, shown as the thin vertical lines) with no baryon bias (b b = 1), black hole to stellar 
mass ratio jbh = 0.0015 and dark matter scaling relations from Sec. [2] The top, solid and dotted, curves represent the total galaxy assuming either Nav04 or 
JS-1.5 dark matter, respectively. The arrows (upper right) indicate the lower limits to the dark matter mass according to Loewenstein & White (1999), given 
our total mass at those radii. 



The figure also indicates that our model produces Nav04 dark 
matter contents that at R = 5 R e are co nsistent with the joi nt kine- 
matics / X-ray modelling of Loewenstein & White 1 1999), while 
our Nav04 dark matter mass frac tion at R = R e is 1 .6 times 
lower than the lower limit found bv lLoewenstein & White! (but our 



predicted JS-1.5 dark matter mass is precisely the lower limit of 
lLoewenstein & White!) . 

Although the gas component begins to dominate the stellar 
component at 10 R e , as expected, the gas does not influence the 
line-of-sight or aperture velocity dispersions: although its influence 
on the line-of-sight velocity dispersion dominates that of the stellar 
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Figure 5. Radial profiles of cumulative dark mass (thin curves, blue in the 
electronic version of the Journal) and stellar line-of-sight velocity disper- 
sion (thick curves, red in the electronic version of the Journal) for the NFW 
(dashed curves) and JS-1.5 (dotted curves) dark matter models, both nor- 
malized to those of the Nav04 dark matter model. Same parameters as in 
Figurel4l 
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Figure 6. Radial profiles of cumulative fractions of gas (long dashes), stars 
(short dashes) and total baryons (solid curve) for a Nav04 dark matter 
model (top), and the cumulative gas/star ratio (bottom). Also shown is the 
baryon fraction with the JS-1.5 model (dotted in top panel). The arrows in- 
dicate the virial radii. Thick and thin curves and arrows are for Tb = 100 
and 390, respectively. Other parameters are the same as in Figurel4l 



component at 60 Re , by that large radius the dark matter component 
is fully dominating the stellar velocity dispersion. 

Given the strong dominance of stars (and the central black 
hole) over the dark matter component, the differences in aperture 
or line-of-sight velocity dispersion profiles between the inner dark 
matter slopes of a — 1 and a — 3/2 are very small, and basically 
indistinguishable from observations. 

The line-of-sight velocity dispersion using the Nav04 profile 
is also indistinguishable from the other two dark matter profiles. 
However, the Nav04 profile produces slightly (5%) lower total ve- 
locity dispersions at R = 5 R e . This is illustrated in Figure [5] 
which plots the mass and line-of-sight velocity dispersion profiles 
normalized to those of the Nav04 model. At radii around r v /8 (cor- 
responding to 10 R e for our standard set of parameters), the NFW 
velocity dispersion is up to 6% larger than that of the Nav04 model, 
which fits much better the density profiles found in dissipation- 
less cosmological A^-body simulations. The corresponding over- 
estimate of the cumulative mass is almost 20% at r v /8 (10 R e ) , 
and the same effect is visible in Figure 2 of iNavarro et"all <2004 
(with the virial radii given in their Table 3). The maximum velocity 
dispersion ratio in Figure|5|indicates that the lower stellar velocity 
dispersions obtained with the Nav04 dark matter model relative to 
its NFW counterpart are not caused by the convergent dark matter 
mass profile of the former at large radii, but by the 20% difference 
in mass profiles at 10 R e . 

The conclusions from Figure|4]are unchanged if we adopt the 
much higher universal total mass-to-light ratio Tb = Tb = 390, 
instead of 100, in particular the stellar component dominates the 
dark matter out to only 2 R e (instead of 3 R e ) for the line-of-sight 
velocity dispersion. 

Note that th e ratio of aperture velo city dispersion at i? e /8 (the 
aperture used bv ljorgensen et alll 19961) to circular velocity at the 



virial radius, a a , p (R e /10)/V v , is 0.64 when Tb = 390, but as high 
as 1.00 when Tb = 100. Indeed, since the stars dominate the aper- 
ture velocity dispersion measurement, decreasing the total mass-to- 
light ratio decreases V v without affecting much a ap (i? e /8). 

Figure[6]displays radial profiles of the baryon fraction and gas- 
to-star ratio. For our standard mass-to-light ratio at the virial radius, 
Tb = 100, the decrease in star fraction with radius is only partially 
compensated by the increase in gas fraction, so that the total baryon 
fraction decreases with radius for r < r v : for the Nav04 model, ft 
decreases from unity at small radii to 0.87 at R e , 0.76 at 2 R e , 0.53 
at 5 R e , and down to a minimum of 0. 132 at R — 68 R e — 0.95 r„, 
and then slowly rises by 0.1% from there to the virial radius. These 
trends are very similar for the NFW dark matter model, while for 
the JS-1 .5 dark matter model, the baryon fraction actually increases 
from 0.87 at i? e /100 to 0.90 at R e /13, then decreases to 0.80 at 
R e , 0.70 at 2 7i e , 0.49 at 5 R e - For the universal mass-to-light ratio, 
Tb = 390, the trends are similar, except that the baryonic fraction 
reaches a minimum of 0.10 (for all dark matter models) around 
50 Re, i.e, r v /2, and rises at larger radii. 

The bottom plot of Figure|6|indicates that the gas-to-star ratio 
rises, roughly as r 1,2 from ru/10 000 to r v . If the mass-to-light 
ratio at the virial radius is increased, while the baryonic fraction at 
the virial radius and the total stellar mass are kept constant, then 
the gas component must become relatively more important, which 
explains why the Tb = 390 gas-to-star ratio is larger than the 
Tb = 100 gas-to-star ratio, by a factor 390/100 = 3.9. 

4.3 Can one weigh the dark matter component? 

As we have seen in paper I and Sect. 14.21 the dark matter contributes 
little to the inner regions of elliptical galaxies, and hence we need to 
focus on the outer regions to be able to weigh the dark component. 
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Figure 7. Line-of-sight stellar velocity dispersions of our 4-component el- 
liptical galaxy (eq. IA9I 1 for our slightly anisotropic (eq. 1601 with r a = 
0.18 r„) NFW (dashed), JS-1.5 (dotted) and Nav04 (solid) dark matter 
models, with total mass-to-light ratio within the virial radius T g = 100 
and 390, increasing upwards (respectively black and red in the electronic 
version of the Journal). Same other parameters as FigureRI 



Figure illustrates the effect of the mass of the dark matter 
component on the line-of-sight velocity dispersion profiles. While 
the steep inner slope of the density profile of the JS-1.5 model al- 
lows it to have a non-negligible effect on the line-of-sight veloc- 
ity dispersions at small projected radii, the opposite is true for the 
NFW and Nav04 models: the inner line-of-sight velocity disper- 
sions are completely independent of the total mass of the galaxy. 
Moreover, at large radii with Tb = 100, it is difficult to distin- 
guish between the three dark matter models. Note that relative to 
the Nav04 dark matter model, the NFW model has no effect on 
the line-of-sight velocity dispersions measured from 0.01 to 0.1 R e 
(even when no central black hole is present), while the JS-1.5 com- 
ponent has a small effect, increasing with mass. 

Figure[S|shows again the effect of the mass of the dark mat- 
ter compon ent, this time allowing the v ery low mass-to-light ratio 
deduced by Romanowskv et al. 1 2003 ) (see Sect. 13. it . which re- 
quires bt > 1 (eq. 1291 ). The figure shows that, at R/R e = 2, 
3, 4, 5, and 6, the velocity dispersion increases by 26, 33, 39, 43, 
and 46kms~ 1 when the mass-to-light ratio at the virial radius is 
increased from Tb = 33 to Tb = 390. In other words, a line- 
of-sight velocity dispersion measurement at 4R e with 20kms _1 
measurement error implies an uncertainty on the total mass-to-light 
ratio of a factor greater than 3\ At higher radii, less precision is re- 
quired on the velocity dispersions, but these are more difficult to 
estimate as they are smaller. These conclusions are the same if we 
adopt purely isotropic models. This three-fold uncertainty becomes 
even larger if we do not assume a precise form for the velocity 
anisotropy profile. 

Also shown in Figure[8|is the case of no dark matter (for which 

= /bTg d = 14). The dispersions now fall to quite low val- 
ues. However, this 'no dark matter' case is extreme, as its baryonic 
fraction is unity. 
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Figure 8. Same as Figure^with Nav04 dark matter, Tb = 33, 100, and 
390, increasing upwards (blue, black, and red in the electronic version of 
the Journal), as well as no dark matter (Tb = /(, T^ d = 14) case (lower 
curve, green in the electronic version of the Journal), for which Tb = 
fb T^ d = 14. The baryon bias is b b = 1 (solid curves) and b b = 1.5 
(dashed curves). The case b\> = 1 is not shown for Tb = 33 because the 
gas fraction would then be negative (see eq. 1281 ). Same other parameters 
as in Figure |4| 



How much do these conclusions depend on the galaxy lumi- 
nosity? Figure|5|displays the dispersion profiles normalized to the 
values at 0.1 i? e (at which radius the dispersion profile is fully dom- 
inated by the stellar component, as seen in Fig. |4j- The velocity 
dispersion profiles have a strikingly similar shape, with relative dif- 
ferences less than 6.5% for R < 5R e , and interestingly, less than 
3.5% for R ~ 2.5 R e . Also, the slope at R > i? e varies with lumi- 
nosity. 

Reverting now to our standard luminosity Lb = Xb — 
1.88 x 1O 1O L0, we look in more detail at how the velocity dis- 
persion at a fixed number of effective radii scales with the total 
mass-to-light ratio. Figure [Tol shows that the velocity dispersion at 
2 and 5 R e rises slowly with mass-to-light ratio, as a power-law of 
slope 1/13 and 1/6, respectively. In other words, the virial mass- 
to-light ratio increases very sharply with the measured line-of-sight 
velocity dispersion at a fixed number of effective radii. Also, at 
fixed measured velocity dispersion, the standard isotropic NFW 
model produces significantly smaller mass-to-light ratios than the 
anisotropic Nav04 model. 

This is quantified in Figure EH At 5 R e , going from the 
isotropic NFW model to the isotropic Nav04 model, the inferred 
mass-to-light ratio is ~ 1.4 times higher (left 'Nav04-iso' curve), 
and this factor is roughly independent of the measured velocity dis- 
persion. Indeed, as seen in Figure [5] at 5 — 20 R e , the velocity 
dispersion of the NFW model is a few percent larger than that of 
the Nav04 model. One therefore needs a higher mass-to-light ratio 
to reproduce a fixed value of a\ OB with the Nav04 model in com- 
parison with the less accurate NFW model. Moreover, going from 
the isotropic NFW (Nav04) model to the slightly anisotropic NFW 
(Nav04) model of equation 1601 . the inferred mass-to-light ratio 
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Figure 9. Normalized line-of-sight stellar velocity dispersion for our 
slightly anisotropic (eq. 1601 with r a = 0.18 r v ) model with Nav04 dark 
matter, with T s = 100 and L B = 3.2 X 10 9 (dotted cun<e), 10 10 (solid 
curve), and 3.2 X 10 10 Lq (dashed curve), with Sersic shape m and effec- 
tive radius R e taken from equations 1161 and 1151 - respectively, and other 
parameters as in Figurel4l 
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Figure 10. Line-of-sight stellar velocity dispersion of our 4-component el- 
liptical galaxy at 2 (upper curves) and 5 (lower curves) effective radii, as a 
function of total mass-to-light ratio at the virial radius (or its bias relative to 
the universal value, upper axis). The dotted, short-dashed, solid and long- 
dashed curves represent the isotropic NFW, isotropic Nav04 and anisotropic 
(eq. 1601 ) Nav04 models, respectively, the latter two with r a /r v = 0.18 
and 0.018, respectively. Same parameters as in Figure in particular no 
baryon bias (6^ = 1), which limits the virial mass-to-light ratios as in equa- 
tion l28l . 



Figure 11. Ratio of the mass-to-light ratio inferred from a general model to 
the mass-to-light ratio inferred from the isotropic NFW model, as a function 
of the measured line-of-sight velocity dispersion at 5 (left) and 2 (right) ef- 
fective radii for a continuous range of virial mass-to-light ratios. The short- 
dashed, solid and long-dashed curves represent the ratios to isotropic NFW 
mass-to-light ratios for the isotropic Nav04 and anisotropic (eq. 1601 ) Nav04 
models, respectively, the latter two with r a /r v = 0.18 and 0.018, respec- 
tively, while the dotted cun'e shows the analogous ratio for the anisotropic 
NFW model (eq. 1601 with r a = 0.18 rv). Same parameters as in Figurel4l 



is 1.1 to 1.2 times larger (increasing for smaller measured veloc- 
ity dispersions), again as expected since radial anisotropy causes 
lower values of <ti os at large radii (see Fig.[3j. The combined ef- 
fect of dark matter model and velocity anisotropy is displayed in 
the upper curve of Fieure fTTI and indicates that line-of-sight ve- 
locity dispersion measurements lower than 135 kms -1 for an L* 
elliptical galaxy imply an underestimation of the mass-to-light ra- 
tio within the virial radius of a factor of roughly 1.6, when using 
the isotropic NFW model instead of the best slightly anisotropic 
models arising from dissipationless cosmological N-body simula- 
tions. Similarly, the velocity dispersion at 5 R e produced by the 
Nav04 dark matter model with the strong radial stellar anisotropy 
(r a = 0.018 r v ) seen by \Dekel et al\ in merger simulations is the 
same as that produced by an isotropic NFW dark matter model with 
a virial mass-to-light ratio 2.3 times lower. 

Similar patterns are found at 2 R e but with considerably 
smaller effects. Note that a positive baryon bias (bt > 1) will pro- 
duce lower velocity dispersions at a given radius for a given mass- 
to-light ratio (Fig.[8j, hence the inferred mass-to-light ratio will be 
increased relative to the value obtained for 65 = 1. 

How far out should one measure velocity dispersions? Ta- 
ble Q compares the uncertainties at 2 and 5i? e , from the two ef- 
fects we discussed above: 1) the observational uncertainty on the 
velocity dispersion as inferred from the slopes of the curves in Fig- 
ure^] and 2) the uncertainty caused by uncertain total (or equiv- 
alently dark) mass and anisotropy profiles, as inferred from Fig- 
ure \H\ where we assume the velocity dispersion profile given in 
Fieure flOl and guess the measurement errors (see also Sect.|5] be- 
low). For assumed measurement errors on a\ os of 5 km s _1 at 2 R e 
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Table 1. Relative uncertainties on the mass-to-light ratio at the virial radius 



R a v Sa, (") ( S -l) ( 5 -l) 

V /obs \ T / mod V T / t. 

(kms -1 ) (kms -1 ) 
2R e 130 5 0~9 L0 lX~ 

5 R e 106 20 1.8 1.4 2.3 

5 R e 106 10 0.7 1.4 1.6 



ity of the galaxy also affects the slope of the line-of-sight velocity 
dispersion profile. 

The larg e range of mass-to-l ight ratios at the virial radius that 
we infer from Roman owskv et all 's 15 orbit solutions (Fig.0 sug- 
gests that the detailed orbit modelling provides little constraint on 
the dark matter content of ellipticals within their virial radii. 



5 SUMMARY AND DISCUSSION 
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Figure 12. Logarithmic slope of the radial profile of the line-of-sight stellar 
velocity dispersion a a function of the total mass-to-light ratio at the virial 
radius. Same line types as in Fig.[l0| 



and 20 km s" 1 at 5 R e , we end up with a total factor of 2.3 at 2 R e 
and of 3.3 at 5 R e , the latter factor being reduced to 2.6 if the veloc- 
ity dispersion at 5 R e can be measured with 10 km s _1 accuracy. 

Of course, one can do better by combining the measurements 
of crios at various radii between 1 and 5i? e - A quantitative as- 
sessment of combining measurements is beyond the scope of the 
present paper. However, inspection of the lower left panel of Fig- 
ured indicates that the logarithmic slope of the line-of-sight veloc- 
ity dispersion profile is roughly independent of the limiting radii, 
when these lie in the interval [1, 20] R e - The middle right panel of 
Figure |4] shows that this is not the case for the (cumulative) mass- 
to-light ratio: the slope rises with radius. In view of the possible 
1.3 (or more) underestimation of the effective radius of NGC 3379, 
any error on estimating the effective radius of a galaxy will lead to 
biases in analy ses based upon the mass -to-light ratio gradient, such 
as the study of N apolitano et al J 12005). 

It is much safer to base one's analysis on the velocity disper- 
sion gradient. Fieure IT2l shows that the logarithmic slope of the 
line-of-sight velocity dispersion profile decreases linearly with the 
log of the mass-to-light ratio at the virial radius. Alas, the direct 
interpretation of this logarithmic slope is complicated by the un- 
certainties in the dark mass and stellar anisotropy profiles: for a 
given slope of the velocity dispersion profile, these uncertainties 
lead to a factor 3 uncertainty on the mass (or mass-to-light ratio) 
at the virial radius. Moreover, as noted earlier (Fig.[5) the luminos- 



In this paper, a 4-component model of elliptical galaxies has been 
built to compare the predictions of cosmological TV-body simula- 
tions with the observations of elliptical galaxies. The inner regions 
are dominated by the stellar component, unless the inner dark mat- 
ter density profile is as steep as r _3//2 , which seem s to disagree 
with t h e latest cosmologic a l TV-body simulation s ofJStoehr^yilJ 
1 2002 ). IFower et all 12003]) . INavarro et alJ 120041) . iDiemand et al J 
(2004a), and Stoeht (2005). Therefore, there is little hope in con- 
straining, through the analysis of internal kinematics, the inner 
slope of the dark matter density profile unless it is as large as 3/2 
in absolute value. 

The dark matter component is found to become important 
at typically 3-5 effective radii. At these radii, the galaxy surface 
brightness is low and one requires the sensitivity of an 8m-class 
telescope, and the required precision on velocity dispersion of say 
5 kms -1 imposes a decent spectral resolution (with R = 1000, 
a signal-to-noise ratio S/N > 10 should lead to lOkms^ 1 preci- 
sion on observed velocity dispersions — L. Campbell, private com- 
munication). For example, according to the Exposure Time Cal- 
culator at the ESO Web site 7 , observations of the giant elliptical 
galaxy NGC 3379, using the the Sersic surface brightness profile 
we fit (see Sect. 13. U . on the VLT-UT1 (Antu) 8m telescope us- 
ing the FORS2 multi-slit spectrograph at R = 2000, one can 
reach S/N = 10 (per spectral pixel and per arcsec) in 3 hours 
at R = 2R e . Among the 16 slits, one can dedicate say 10 for 
R = 3 R e and the combined spectrum should produce S/N > 10 
with a precision on velocity dispersion of order 5 km s~ . This 
measurement at 3R e should constrain the total mass-to-light ra- 
tio within the virial radius, but to fairly low precision: at 3 R e , one 
has d In Tb /d In a v ~ 9, yielding an accuracy of a factor of 1.5 
only. 

Three physical mechanisms are found to decrease the ob- 
served line-of-sight velocity dispersion at 2-5 effective radii: 

(i) the lower cumu lative mass of the new dark matter models 
such as that of INavarro et alJ 120041) at these radii and slightly 
above; 

(ii) radial velocity anisotropy; 

(iii) a possible excess of the baryon fraction relative to the uni- 
versal value. 

In turn, if we model the observed velocity dispersions at 5 ef- 
fective radii, the combined effect of the new dark matter models and 
velocity anisotropy is to increase the inferred mass-to-light ratio at 
the virial radius by 60% for slight anisotropy, as observed for the 
particles in structures within dissipationless cosmological TV-body 
simulations, and b y a factor 2.4 for the strong anisotropy found by 
lDekeletak. (2005) in simulations of merging galaxies. 

The effect of the velocity anisotropy on the mass at the virial 
radius is a direct illustration of the mass / anisotropy degeneracy. 



http://www.eso.org/observing/etc 
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On the other hand, the analysis of the galaxy kinematics at some 
radius R requires the knowledge of the mass profile beyond R in 
theory, but just a little beyond R in practice, since the mass profile 
entering the expression (eq. IA15I ) for the line of sight velocity 
dispersion is weighted by the luminosity density profile, which falls 
relatively fast. Therefore, the extrapolation of the mass profile out 
to the virial radius depends on the details of the mass profile. 

In their analysis of halos from dissi pationless cosmologi- 
cal iV-body simulations, Sanchis, Lokas, & Mamon (2004) showed 
that the mass profile was recovered to high accuracy beyond 
0.03 r v . The difference here, is that the line of sight velocity disper- 
sion data in ellipticals extends only out to 2 to 5 effective radii, i.e. 
to less than 6% of the virial radius, while in clusters of galaxies, the 
analogous dispersions extend all the way out to the virial radius. 

One may be tempted to explain the very lo w (Tb = 33) mass- 
to-ligh t ratio at the virial radius reported by IRomanowskv et all 
(2003) (which, at face value requires the baryonic fraction within 
the virial radius to be larger than the universal value) by this factor 
of 2.4, which would raise the mass-to-light ratio at r v to Tb w 80, 
and to roughly 100 when we use rio2 for the virial radius instead of 
r2oo- However, this comparison is not fair, because the best fitting 
orbital solutions of IRomanowskv et alj are also anisotropic: their 
orbits are nearly as radial as the stellar orbits at a few R e found by 
iDekel etHll cA. Romanowsky, private communication). 

Furthermore, the low dispersions of pl anetary nebulae ve- 
locities around nearby ellipticals measured bv lRomanowskv et all 
« 80kms -1 , cannot be reproduced with our models: as seen in 
Fig. El we consistently produce velocity dispersions at 5 R e above 
cios = 95kms _1 for Tb > 80. Therefore, resorting to radial 
anisotropy and to the Nav04 dark matter model appear to be in- 
sufficient to explain the low velocity dispersions observed at large 
radii by IRomanowskv et all And yet, we found that some of the 
orbit solutions found by these authors lead to high enough mass-to- 
light ratios (Fig.0, as does a new solution (Romanowsky, private 
communication, and right-most pair of points in Fig, [p. In fact, 
the 2nd most massive among the solutions of Romanows ky et all 
(third pair of points from the right in Fig.0 has Ts(rio2) = 70, 
and if a Nav04 dark matter model were used instead of the NFW 
model, one would get roughly 40% larger masses at the virial radius 
(see Fig.EJ, i.e. 100, as rou ghly expected from cosmological mod- 
elling iMarinoni & Hudsonll2002l:lYang et alj|2003l) . Note that the 
more massive solutions of Romanows ky et alJ have low dark mat- 
ter concentration pa rameters, co mpared with the expect ations of 
Bullo ck et alJ fcOOll) . revised bv lNapolitano et all feOOSl) . but one 
solution (T_b.io2 = 70, c = 10) is consistent with the predictions 
from dissipationless cosmological simulations. 

Dekel et al. suggest other ways to bring down the velocity dis- 
persions: a steeper density profile for the planetary nebulae mea- 
sured by IRomanowskv et all and viewing triaxial galaxies along 
their minor axis. Another possibility is an incomplete dynami- 
cal equilibrium of elliptical galaxies, caused by small residual 
time-variati ons, aff ecting both the orbit modelling performed by 
IRomanowskv et ail and the Jeans kinematical modelling, as pre- 
sented here. Detailed analyses of kinematical observations and of 
merger simulations are in progress to help clarify the dark matter 
content of elliptical galaxies. 
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APPENDIX A: LINE-OF-SIGHT VELOCITY DISPERSIONS FOR ANISOTROPIC MODELS 
Al Radial velocity dispersion 

The general solution to the Jeans equation is 

^)^W = y^) / f(s)e(s)^^lds, (ad 

where / is the solution to 

^=2/?(r), (A2) 
dmr 

i.e., 

ds 



f(r) = /(n) exp 



20(a) 



s 



(see also van derMarejl994h . 

For isotropic orbits, the Jeans equation trivially leads (/ = 1) to 

For radial orbits, one has iLokas & Mamorl2003l) 

„ poo 

£(r) cr 2 (r) = -± / fMds. (A4) 
Also, with the Osipkov-Merritt anisotropy (eq. 1611 ) one gets with / = r 2 + v\ 

^)°l(r) = -jiLj d S - (A5) 

Similarly, with the anisotropy profile of equation 1601 . shown as the solid curve in Figure |2| one gets with equations JA1> and <A2t and 

/ = r + r a 

£(r) a 2 r (r) = [°° I M (^±pi) ds . (A6) 

r + r a J r \ s 2 J 

A2 Line-of -sight velocity dispersion 

A2.1 General formulas 

Projecting the velocity ellipsoid along the line-of-sight, one finds that the line-of-sight velocity dispersion is (Binnev & Mamon 1982) 

i -i( R ) aL(R) = r - i? 2 r iifi^ . ( A7) 

Inserting equations <AH and <A2t into equation <A7t and inverting the order of integration, one obtains the general expression 
HR)"UR) _ 2 r f*M [ s 1 rdr 2 [°° f£M H df/dr dr 

g - 2 J R — ds J R 77^W- R 1 R — ds J R —V^W (A8) 

The inner integrals of equation <A8> are closed for only several realistic 8 forms of /: 

(i) / = r 2 :/3 , whichimpliesP = est; 

(ii) / = (r + a) a , with a = 1 or 2, which implies (3 = (a/2) r/(r + a); 

(iii) / = (r 2 + a 2 )", aga in with a = 1 or 2, which implies /3 = ar 2 /(r 2 + a 2 ) — hence, a = 1 corresponds to Osipkov-Merritt 
( Osip kovi 1 979UMerrittl 1985h anisotropy. 

For the anisotropy profile of equation 1601 . one obtains from equation <A8> . for R ^ r a , after inverting the order of integration, with 
plenty of algebra 



I(R)aL(R) f°° ( r + r a \ .-i/r\... „, . . , R 2 /~°° Vr 2 - R 2 

g = 2 J R J cosh {R) £ ^ M ^ dr+ VT^W J R ~ P~ 

f r + r a \ c _j / r a r/R + R \ 
\ r 2 ) \ r + r a J 

8 we reject / = \Jr 2 — a 2 and / = 1/V r 2 — « 2 , as both lead to unphysical /3 > 1 near r = a. 



-sgn^-lj r a ° / (— "K i " iMD.l/indr. UV» 
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where C 1 (x) — cosh 1 x for r a > R and C x (x) — cos 1 x forr„ < R. 9 For R = r a , we similarly obtain 

I{ra)a? os (r a ) 



G 



(^) cosh "(£) *(r)M(r)dr-±|" (^^) dr . 



(A10) 



In the limit of isotropy (r a — > oo), we recover (after some algebra on eq. IA9I ) the standard line-of-sight velocity dispersion for isotropic 
systems (Prugniel & Simien 1997): 



I(R)aL(R)=2G 



\/r 2 - R? 

2 £(r)M{r)dr 



(All) 



which can also be easily obtained through equations <A3> and <A7> . again after inversion of the order of integration. 
For Osipkov-Merritt anisotropy (eq. 1611 ) one finds, in a similar fashion, 



I(R)aL(R) _ 2r 2 a + R 2 

(r 2 + R 2 ) 3 / 2 JI; 



r 2 +r 2 \ i lr 2 — R 2 

' tan" 1 \ I - „„ i{r) M(r) dr 



R 



G (r 2 + R 2 ) 3 ' 2 L V r 2 / V d+ff 

For constant anisotropy orbits, the same procedure yields 



r 2 +R 2 



Vr 2 -R 2 



£(r) M(r) dr 



(A 12) 



R 



1-2/3 



r(/ r(/3)' /2 ' ) ^ 'lr).U(/-)dr 4- I I BGi 2 /r..l- -1./2. J./2)/(r) . I/O), Ir 



B(R 2 /r 2 ,f3- 1/2,1/2) l{r)M{r)Ar 



(A13) 



where a, b) = t a 1 (1 — t) 6 Mt is the incomplete beta function. For purely radial orbits, one can obtain the simpler formula 



G 



2R 



£(r) M(r) dr - 



Vr 2 - R 2 a , ^ , J 1 
e(r) M{r) dr- — 



IR JR " JR 

Therefore, for our choice of simple anisotropy profiles, we can generally write 
I(R) aL(R) = 2 G J™ K (1 t(r) M(r) * , 
where the kernel _K" can be expressed for our five anisotropy models as 



(f) t(r)M{ 



r) dr . 



1 2/3-1 

2 U 



7T 1 / 1 U , - 

■ u — — \ 1 ?r — — sin 



if(u,ti a )=< «4 + l/2 /ll 2 +<\ _! /-U 2 -l 

' tan 



(^i + i) 3/2 V u 



1/2 



Ut + l U 2 a + 



jVi-iA 



JZ^ /^ + ( 1 + M cosh - lu 
ug — 1 v u V w / 

/ -n — 1/2 / u a \ 1 (u a u + l 



(cst-/3) , 



(isotropic) , 



(radial) . 



(Osipkov-Merritt) 



(eq. GJB,u / 1) , 



1 + -1 cosh 1 u - - I 7 



1/8 
6 \u 



it- 1 



it / b \u / y u + l 

where _B is the incomplete beta function, C _1 (x) = cosh -1 a; for u a > landC _1 (x) = cos -1 x for m„ < 1 



(eq. |60|, u a = 1) , 



(A14) 



(A15) 



(A 16) 



A2.2 Application to cold dark matter profiles and a Sersic luminosity profile 

The total mass at a given radius is the sum of masses from the dark, stellar, gas and black hole components: 

M(r) = M t M{r/a s ) 



(A17) 



9 Here and below, we use the positive definitions of cos 1 x and cosh 1 x. 
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M(x) = U M d (r]x) + /. + f g Z 3( ^ X) + .gBH /. , (A18) 

L 3 (x v ) M g (£x v ) 

where we used equations i 10t . 1 1 81 . J32I . J37> . 1381 . <47> . and 1481 - as well as the component fractions at the virial radius, fi, from equations 
1231 - 1261 . where x — r/as, x v = r v /as, t] = as /ad, £, = as/r c = q/b m , the deprojected luminosity L3 of the Sersic profile is given in 
equation II li . and the dimensionless gas profile M g is given in equation 1461 for r ^ r v or I54i for r > r v . One can set /;,=/» = T* = 
if one wishes to study the case of an NFW, JS— 1.5, or Nav04 potential, i.e. assuming that the total matter density profile is the NFW, JS— 1.5, 
or Nav04 density profile. 

Using equations 0, ^6), Q, {8}, IA15> . IA17> and IA18> the line-of-sight velocity dispersion of component i can be written 

°L,i( R ) = rwi 2m ^ 1 — ^ ex P ( xVm ) * I?) fc) ^) — > < A19 > 
l[(6—p)m\ as y Jx ^ X X ' x 

where X = R/as, x a = r a /as, and the kernel K is given in equation IA16> . Then of os (7?) = J^. of OSi i(-R). 
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